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1. Introduction 

Recent numerical work provides evidence that Wilson loops in SU (N) gauge theory in two, 
three and four dimensions exhibit an infinite -N phase transition: Dilated from a small size to a large 
one, the eigenvalue distribution of the untraced Wilson loop unitary matrix expands from a small 
arc to the entire unit circle [1, 2]. This transition, which was discovered by Durhuus and Olesen 
[3], is in the universality class of a random multiplicative ensemble of unitary matrices [4]. In 
the following, we will relax the unitarity constraint and focus on a multiplicative random complex 
matrix model introduced and solved in [5]. 

2. Basic random complex matrix model 

We define a sequence of n independent and identically distributed N xN matrices 

Mj = exp(eCj) , j = l,...,n (2.1) 
with normalized Gaussian probability distribution 

P(C)dli(C)=e- mrC * c n " MQj) (2-2) 

1<U<N K s " ' 

dRe{Cij)dlm{Cij) 

which is invariant individually under Cj — ► C*,—Cj and Cj — ► UjCjUj, Uj £ U(N). We do not 
restrict the trace of Cj here since it turns out that requiring detMy = 1 has no effect on the saddle- 
point analysis in the large-A^ limit. We are interested in the distribution of the product 

n 

W n = ]\Mj (2.3) 

7=1 

in the limit n — > °°, e — ► with the parameter t = e 2 n held fixed, t > can be interpreted as a 
diffusion time. This model is almost identical to that of [5]. 

To derive a closed formula for the entire distribution of W for general ,/V is difficult; however, 
similarly to the unitary model, partial information about the distribution of eigenvalues can be 
obtained from the averages of characteristic polynomials. In the following, averages over all Cj are 
denoted by (...). (det(z — W n )) carries no information, since 

(dct(z-W n )) = (z-l) N . (2.4) 

The first non-trivial case is 

Q(z,z*) = {\det(z-W n )\ 2 ). (2.5) 

Applying large-Af factorization to Q (i.e., assuming that the average of the product can be replaced 
by the product of the averages) would result in holomorphic factorization, 

Q(z,z*) - (det(z-wg>(det(z-wg*> = \z-l\ 2N , (2.6) 

and all eigenvalues seem to have to be unity. This factorization is expected to hold only in two 
distinct regions: inside a circle around zero with radius p < 1 and outside a circle of radius p _1 > 1. 
Therefore, the surface eigenvalue density is restricted to the annulus p < |z| < p 1 . The aim of the 
following is to calculate Q as a function of t. We shall find that the domain of eigenvalues becomes 
multiply connected at a critical t = t c , in agreement with [5]. 
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3. Saddle-point analysis 



The first step in the procedure is to disentangle the non-abelian product of matrices defining 
W n . To this end, we introduce 2nN pairs of Grassmann variables { xj/j, Yj,XjiXj}j=\ „ an d find 

|det( z -W„)| 2 = J fl dyjdVjdXidXj^ : '" V¥ - e ° x x e r ■ " " " x X ' x ' (3.1) 

with z = e na and the convention that Yn+i = Vi> etc - 

Now, the integrals over the matrices Cj factorize and can be done explicitly to sufficient accu- 
racy in £. The following equalities ought to be understood in the sense that they hold up to terms 
which vanish as n — > <*>, e — > at t = e 2 n fixed. An expansion of M = exp(eC) to linear order in e 
is sufficient because the next term does not contribute, 



(e 



-W-xx'-jiVx'xV 



(3.2) 



Introducing scalar complex bosonic multipliers ^, j = 1, . . . ,n allows for a separation of the 
quartic Grassmann terms into bilinears, 



n J 



(3.3) 



where the integration measure is djj,(Q = dRe £ dim £. After shifting some indices of Grassmann 
variables, this leads to 



Q(z,z*) 



N Y f f\[dWjdWjdXjdXjdli(Cj)}e- Nr J=^ 
nj J j= i 



(3.4) 



As a result, averages over complex matrices Cj are replaced by averages over complex numbers £y. 
Carrying out the integrals over the Grassmann variables makes the dependence on N explicit, and 
we are left with 



where 



A = e f7 r t -1, D=A t , B = -eZ, C = £Z j 



(3.5) 



(3.6) 



with n x n matrices 



/0 1 
1 


\10 



• o o\ 

• 

• 1 

• 0/ 



and Z = diag(Ci,C2, •••,£») ■ 



(3.7) 
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The large-./V limit leads to saddle-point equations trivially satisfied at £j = for all j = 1 , . . . , n 
since the Q enter only bilinearly in 



det 



A B 
CD 



|detA| 2 det(l +e 2 Z t A~ 1 Z(A t ^ 1 ' 



Where this saddle dominates we obtain 



\2N 



(3.8) 



(3.9) 



and W„ could be replaced by a unit matrix, which means that there are no eigenvalues at any z ^ 1 
in the complex plane. 

Comparison with numerical simulation shows that the trivial saddle point is always dominating 
whenever it is locally stable. At the boundary of the local domain of stability one has a transition to 
regions with non-zero surface eigenvalue density. To determine this boundary we need to expand 
the integrand around = 0, 



det 



A B" 



|det(i- e ff r+)|V 2 ^=iW< 



(3.10) 



with Kji = \(A )ji\ ■ The matrix T implements cyclical one-step shifts: T" = 1. Hence, 



? CT rt - 1 1 



— V e- SG T s 

-no L-i 



s=\ 



implying that K is a circulant matrix. Its eigenvalues are found to be 

2 



1 



7=1 



1 - e- na 
The condition of local stability, 



1 



1 - e- na 



1 „-«(ff+C7*) . 



Re(-1 + £ 2 A^) <0, k= l,...,n 



(3.11) 



(3-12) 



(3.13) 



is strongest for k = n, and consequently the region of local stability of the trivial saddle point is 



1 



1 



k-i 



< i 



i 



(3.14) 



Taking the limit (n — > °°, £ — > with t = e 2 n kept fixed) gives 



2\z— 1| 2 loglzl 



(3.15) 



in agreement with [5]. 

In contrast to (3.14), the last inequality is invariant under z — ► The transition occurs when 
the inversion invariant point z = — 1 on the unit circle first enters the domain of eigenvalues. The 
condition for vanishing eigenvalue density at a point z on the unit circle reads 



?<|z-l| 2 for |z| = l. 



(3.16) 
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With z = e l V, this is equivalent to 

cosi//-<l-^. (3.17) 

When t > 4 this condition is clearly violated for any Consequently, the region of non- vanishing 
eigenvalue density contains the whole unit circle. On the other hand, for t < 4 (3.17) will be fulfilled 
for some points on the unit circle. We see that the unit circle contains an arc, centered at z = — 1 
with endpoints at angles y satisfying cos(y) = 1 — 4, which lies completely in the domain of zero 
eigenvalue density. The domain of non-vanishing eigenvalue density becomes multiply connected 
for t > 4. 



4. Generalized model 

The basic random complex matrix model can be generalized to interpolate between the cases 
where the individual factors are unitary or hermitian. To this end, we write each matrix C as a 
linear combination 

C = Hi+iH 2 with h\ 2 =Hi.2 (4.1) 
and introduce two weight factors (0\ .2 > in its probability distribution, 

P{C) = J y e A^ + ^). (4.2) 

Setting (0\ = CO2 = \ we get back to the basic model, and for (0\ — > we are in the unitary case, 
where the spectrum is restricted to the unit circle. 
For (0\ ^ G>2, 

(det(z-W n ))=J(z) (4.3) 

is no longer equal to (z — 1)^, but to a more complicated polynomial in z. The polynomial J(z) 
is completely determined by the two-point function of the matrix C, which depends only on the 
difference o>i — (0\ . Therefore we can simply go to the unitary case with (0\ = 0, for which the 
polynomial J(z) is known from previous work on products of unitary matrices [2], and absorb the 
dependence on (O2 — (0\ in t by a rescaling, 

t± = *(G>2±£»i), f+>|/_|. (4.4) 
Again, where holomorphic factorization, 

!2(z,z*) = (|det(z-^)| 2 ) = |/(z)| 2 , (4.5) 

holds we have no finite surface charge density. To determine for which values of z this formula 
no longer gives the correct answer we apply a strategy similar to the one presented in Sec. 3. The 
only complication is that more quartic Grassmann terms need to be decoupled. Thus, integrals over 
real multipliers r\j, j = 1, . . . ,n have to be introduced in addition to the complex noise factors 
£j. However, the complex variables again enter only as bilinears ^jQ, leading to a trivial saddle 
point £j = at large N. The eigenvalue density vanishes where this saddle dominates because the 
remaining E, and 77 integrals factorize, resulting in holomorphic factorization for Q. Therefore, 
only the structure of the £ saddle determines if z is in a chargeless region, but the dominating 
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saddle points of the E, , v\ integrals affect the local stability of the trivial saddle point at £j = 0. The 
condition for local stability and vanishing eigenvalue density in the desired limit (e — ► 0, n — > oo 
with t = e 2 n kept fixed) is eventually found to be 

l> J + iJf^! - ( 4 - 6 ) 
2\z— lp log |z| 

This inequality is similar to Eq. (3. 15) for the basic model, with t replaced by t + and z replaced by 
z, where the variable £ is related to the original variable z via 



ze-'- 1 ' with u = - [j—[J ■ ( 4 - 7 ) 

The topological transition in the domain of non-vanishing eigenvalue density thus occurs at t+ = 4 
for the generalized model. 



5. Numerical results 

Since we did not explicitly identify the competing non-trivial saddle points, nor determine 
the global stability of the trivial saddle, to establish the transition more evidence is needed, which 
we provide by numerical simulations. It is more convenient numerically to work with the linear 
model M = 1 + eC introduced in [5], instead of our exponential model M = exp(eC). Repeating 
the analysis for this model, it turns out that the boundaries of non-vanishing eigenvalue density 
are equivalent up to a scaling by a factor of exp(— 1-/2), which is equal to unity for the basic 
random complex matrix model. The following figures show results of numerical simulations for 
the linear model performed with matrix dimension ,/V = 2000 and n = 2000 factors in each matrix 
product for ensembles consisting of about 500 product matrices. Figure 1 corresponds to the basic 
model. The boundaries obtained from the stability analysis, indicated by the solid lines in red, 
are in very good agreement with numerical data given by the blue points. Figure 2 shows results 
of numerical simulations for (0\ = ^ , (02 = \ (data points as well as boundaries are scaled with 
the corresponding factor of exp(— 1-/2)). The topological transition occurs at t = 20/3, which 
corresponds tof + = y Q + -^j = 4 in agreement with the prediction. 
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Figure 2: Same as Fig. 1 but for (0\ = ^, ©2 = \ and t = 5 (left), t = y (middle), t = 10 (right). 



6. Conclusion 

Our main objective in re-analyzing the model of [5] is our conjecture that it is a universal 
representative of the large-jV phase structure of classes of complex matrix Wilson loops. A gen- 
eralization of the probability distribution allows for an interpolation between the cases where the 
individual factors are hermitian and unitary. We confirm, analytically and numerically, that a topo- 
logical transition in the domain of non-vanishing eigenvalue density indeed occurs. Like in the 
unitary case, we would like to study in more detail the approach to infinite N and see what the 
matrix model universal features of this transition are. Since a full analysis keeping the exact N 
dependence is complicated this has not been carried to completion yet. 
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